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Abstract — Computational Anatomy (CA) has introduced the 
idea of anatomical structures being transformed by geodesic 
deformations on groups of diffeomorphisms. Among these ge- 
ometric structures, landmarks and image outlines in CA are 
shown to be singular solutions of a partial differential equa- 
tion that is called the geodesic EPDiff equation. A recently 
discovered momentum map for singular solutions of EPDiff 
yields their canonical Hamiltonian formulation, which in turn 
provides a complete parameterization of the landmarks by their 
canonical positions and momenta. The momentum map provides 
an isomorphism between landmarks (and outlines) for images 
and singular soliton solutions of the EPDiff equation. This 
isomorphism suggests a new dynamical paradigm for CA, as 
well as new data representation. 

I. Introduction 

Computational Anatomy (CA) must measure and analyze 
a range of variations in shape, or appearance, of highly 
deformable structures. Following the pioneering work by 
Bookstein, Grenander and Bajscy [4], [8], [22], the past 
several years have seen an explosion in the use of template 
matching methods in computer vision and medical imaging 
[3], [14], [24], [25], [30], [33], [36]-[40], [51], [53]-[56]. 
These methods have enabled the systematic measurement and 
comparison of anatomical shapes and structures in biomedical 
imagery leading to better understanding of neurodevelopmen- 
tal, neuropsychiatric and neurological disorders in recent years 
[5], [6], [11], [13], [17], [20], [42]-[48], [50], [52], [59]. 
The mathematical theory of Grenander's deformable template 
models, when applied to these problems, involves smooth 
invertible maps (diffeomorphisms), as presented in this context 
in [55], [56], [14], [41], [36] and [33]. In particular, the 
template matching approach involves Riemannian metrics on 
the diffeomorphism group and employs their projections onto 



specific landmark shapes, or image spaces. 

On the other hand, the diffeomorphism group has also 
been the focus of special attention in fluid mechanics. For 
example, Arnold [1] proved that ideal incompressible fluid 
flows correspond to geodesies on the diffeomorphism group, 
with respect to the metric provided by the fluid's kinetic 
energy. In this paper, we shall draw parallels between these 
two endeavors, by showing how the Euler-Poincare theory of 
ideal fluids can be used to develop new perspectives in CA. 
In particular, we discover that CA may be informed by the 
concept of canonical momentum for geodesic flows describing 
the interaction dynamics for singular solitons in shallow water 
called peakons [9]. 

Outline of the paper. Section [H] describes the template 
matching variational problems of computational anatomy, and 
introduces the fundamental EPDiff evolution equation, which 
describes the evolution of the momentum of an anatomy, (a 
collection of landmarks, for example). The singular solutions 
for the EPDiff equation Q with diffeomorphism group G are 
explained in section [HI] They are, in particular, related to the 
landmark matching problem in computer vision. The conse- 
quences of EPDiff for computational anatomy are described 
in section HV1 Conclusions are summarized in section IV1 

II. Variational Formulation of Template 
Matching Problem 

A. Geometrical large deformation setting. 

Grenander [23] pioneered the introduction of group actions 
in image analysis, through the notion of deformable templates. 
Roughly speaking, a deformable template is an "object, or 
exemplar" Iq on which a group Q acts and thereby generates, 
through the orbit I = QIq, a family of new objects. This 
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"template matching" approach has proven to be versatile and 
useful in different settings (image matching, landmark match- 
ing, surface matching and more recently in several extensions 
of metamorphoses [57]) The approach focuses its modeling 
effort on properties of the families of shapes generated by 
the action of the group Q on the deformable templates. Right 
or left invariant geodesic distances on the group Q are the 
natural extension to large deformation of the quadratic cost, 
or effort function, defined for small linearized perturbations of 
the identity element. 

B. Case of non-rigid template matching. 

The optimal solution to a non-rigid template matching 
problem is the shortest, or least expensive, path of continuous 
deformation of one geometric object (template ) into another 
one (target). For this purpose, we have introduced a time- 
indexed deformation process, starting at time t = with the 
template (denoted and reaching the target at time t = 1. 
At a given time t during this process, the current object I t is 
assumed to be the image of the template, I n , through the (left) 
action of a diffeomorphism <fi t : I t = fflo- The attribution of a 
cost to this process is then based on functionals defined on the 
group of diffeomorphisms, following Grenander's principles. 

A simple and natural way to assign a cost to a diffeomorphic 
process indexed by time is based on the following: to measure 
the variation ft+dt — ft, express this difference as a small 
vector of displacement, dtu t , composed with (p t , 

ft+dt - ft = dt u t o (p t 

The cost of this small variation is then expressed as a function 
of u t only, yielding the final expression, 

Cost(t h-> ft) = / £(u t ) dt 
Jo 

with 

— =u t o(p t (1) 
dt 

In the following, the function u t i— » i(ut) is defined as 
a squared functional norm on the infinite dimensional space 
of velocity vectors. This process is a standard construction in 
the Riemannian geometry of Lie groups, in which the con- 
sidered group is equipped with a right invariant Riemannian 
metric. Here, the vector space of right invariant instantaneous 
velocities, u t = (d<j> t /dt) o f^, forms the tangent space at 
the identity of the considered group, and may be regarded as 
its Lie algebra. This vector space will be denoted g in the 
following (as a formal analogy Lie algebra notation). In this 



context, the cost of a time-dependent deformation process, thus 
defined by 

Cost{t^ ft) = f \u t \ 2 B dt, (2) 
Jo 

is the geodesic action of the process for this Riemannian 
metric, and most problems in CA can be formulated as 
finding the deformation path with minimal action, under the 
constraint that it carries the template to the target. We will 
illustrate this below with an important example of a such 
problem, in which the geometric objects are collections of 
points in space (landmarks). Before doing this, we summarize 
in more rigorous terms the process described above, providing 
at the same time the notation to be used in the rest of the 
paper. However, most of the remaining of the paper can be 
understood by referring to the summary paragraph at the end 
of section IlLCl 

C. Rigourous construction. 

Fix an open, bounded subset £1 C R . Following [55], the 
construction is based on the design of the "Lie algebra" g, 
which is in turn used to generate the group elements. (This is 
the converse of the usual consideration of finite dimensional 
Lie groups.) The following construction of g will be assumed. 
Denote by H the set of square integrable vector fields on Q 
with the usual L 2 metric ( • , • ) : H x H — > R. Consider a 
symmetric and coercive operator L : u h Lu G H* = H 
whose domain D(L) contains all smooth (C°°) vector fields 
with compact support in 0. This operator induces an inner 
product on D(L) by (u, w) = (Lu, w). This pre-Hilbert 
space can be completed to form a Hilbert space ( [60]), thereby 
defining g which is continuously embedded in H (Freidrich's 
extension). Suppose, in addition, that g can be embedded into 
C p (fl), the set of p times continuously differentiable vector 
fields on SI, with p > 1. (We call this the p-admissibility 
condition.) Then the following can be shown ( [14], [55]): 
If u t is a time-dependent family of elements of g such that 
J Q \\u t \\ 2 dt < oo, that is, u t G L 2 ([0, l],g), then the flow 
^ = UtO(p t with initial conditions fo(x) = x, x G £1, can be 
integrated over [0, 1], and f\ is a diffeomorphism oftt, which 
is denoted f±. The image of L 2 ([0, l],g) by u i— > (p± forms 
our group of diffeomorphisms, Q, which is therefore specified 
by the operator L. In this setting, equation ([0 has solutions 
over any finite interval, and the infinum should be taken with 
respect to f — ip u a <p for u G L 2 ([0, l],g) such that (pi = 
fi o ip . (Because of right-invariance, it can furthermore be 
assumed that <fro = id.) Moreover, as proved in [14], [55], the 
existence of a minimizer (a geodesic path) is guaranteed. 
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For the inner product (u, w) fl = (Lu, w), the operator L is 
a duality map. In a deformation process 1 1— > t^t such that 

u f is called the (Eulerian) velocity, belonging to g, and Lu t 
is called the momentum, also denoted m t . Note that, because 
L will typically be a differential operator, the momentum Lu 
can (and will) be singular, e.g. a measure, or a generalized 
function. A simple example of such a phenomenon occurs 
in the landmark matching problem described in the next 
subsection. 

This framework for CA is reminiscent of the least action 
principle for continuum motion of fluids with Lagrangian 
t(u) = |||u||g. Note that in the template matching framework, 
I has the specific interpretation of an effort functional for small 
deformations that should be designed according to a given 
application and not follow any existing physical model. How- 
ever, this similarity with ideal fluid dynamics sets the stage for 
"technology transfer" between computational image science 
and fluid dynamics - e.g., Hamiltonian description, momentum 
evolution, classification of equilibria, nonlinear stability, PDE 
analysis, etc. In particular, the least action interpretation of Q, 
is central in the Arnold theory of hydrodynamics [1] and the 
derivation of the geodesic evolution equations falls into the 
Euler-Poincare (EP) theory, which produces the EP motion 
equation [27], [41], 

(J^ + u • v)m+ V*u • m + m(divu) = 0, (3) 

and u = G * m, where G* denotes convolution with the 
Green's kernel G for the operator L. This is the EPDiff equa- 
tion, for "Euler-Poincare equation on the diffeomorphisms". 

Summary: The important consequences of the previous 
construction is that, by measuring the amount of fluid defor- 
mation which is required to morph an object to another, this 
measure being given by (0 where u t is the velocity of the 
fluid deformation at time t and u t |^ = (Lu t ,u t ), L being a 
linear operator, the associated momentum, m t = Lu t satisfies 
the Euler Poincare equation (|3}. This equation is important, 
because it allows to reconstruct the complete evolution of 
the momentum (and hence of the fluid motion) from the 
initial conditions. This property is exploited for the analysis 
of landmark data in [58]. 

D. Landmark matching and measure-based momentum 

The landmark matching problem is an interesting illustration 
of the singularity of the momentum which naturally emerges in 
the computation of geodesies. Given two collections of points 
X\ , . . . , Xk and Y\ , . . . , Yk in Q, the problem consists in 



finding a time-dependent diffeomorphic process (t i— > ip t ) of 
minimal action (or cost, as given by (0) such that tp = id 
and (pi(Xi) = Yi for i = 1,...,K, This problem was 
first addressed in [31], then studied in different forms in 
[7], [10], [18], [19]. Its computational solution relies on the 
key observation that the problem can be expressed uniquely 
in terms of optimizing the landmark trajectories, Q;(i) = 
ift(Xi), for an action, or cost, given by 

s = [ e(Q,Q')dt = l f 'Q'W^QWr'Q'W*, 

J ° J ° (4) 

with notation Q'(t) = dQ/dt. This action, or cost, S = 
J £(u) dt is the time integrated Lagrangian in the least action 
principle, 5S = 0. Its end point conditions are Qi(0) = Xi 
and Qi(l) = Yi, where A(Q) is an Kd x Kd matrix (d is the 
dimension of the underlying space) which may be constructed 
as follows. Let G be the Green's kernel associated to the 
operator L, formally defined by 

v(x)= f G(x,y){Lw){y)dy . 

Let Id be the d-dimensional identity matrix. Then A(Q) is 
a block matrix (Aij(Q),i,j = 1,...,K) with Ay(Q) = 

Denote P(i) = A(Q(t))- 1 Q'(t) = d£/dQ'. Then, the 
optimal diffeomorphism t ip t is given by equation Q with 

K 

u t (x) =^G(x,Q i (t))P i (t). 

i=l 

The corresponding momentum m t is given by the point 
measure 

K 

mt(y) = X)P*(*)%-Qi(*))- 

i=l 

A straightforward extension of this model occurs when the 
landmarks are organized along continuous curves in which 
the indices i, j are replaced by curve parameters (say defined 
over [0, 1]), and 

mt(y) = f P(t,s)8(y-Q(t,s))ds. 

One could also distribute the landmarks along several continu- 
ous curves (the outlines of an image, say). This representation 
of momentum would involve both integrations and sums. 

As we will see, this measure-based momentum is also found 
in other contexts very different from medical imaging. We also 
mention an important variant of this matching problem. This 
variant is a particular case of a general framework, in which 
the diffeomorphic group action is extended to incorporate 
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possible variation in the template itself. The application of 
this theory of "metamorphoses" [57] to the particular case 
of landmark matching simply comes with the addition of a 
new parameter a 2 > and the replacement in the above 
formulas of A(Q) by A° (Q) = A(Q) + a 2 I Nd (I Nd is the 
Nd dimensional identity matrix). This corresponds to geodesic 
interpolating splines introduced in [10]. Note that in this case, 
EPDiff is not satisfied anymore, but may be modified to 
include a non-vanishing right-hand term. 



III. EPDiff and its singular solutions 

A. The EPDiff equation. 

The EPDiff equation is important in fluid dynamics, because 
it encodes the fundamental dynamical properties of energy, 
circulation and potential vorticity. A first example comes with 
choosing the differential operator L as Lu = u, so that m = 
u, and taking incompressible vector fields, so that div u = 0. 
In this setting, equation Q provides the Euler equations for 
the incompressible flow of an ideal fluid [1], [2]. 

Another physically relevant form of the EPDiff equation is 
the evolutionary integral-partial differential equation with 
[26]-[28] 



m = Lu = u — a 2 Au , 



(5) 



so that 



(|u| : 



|Vu| 2 )cfe. 



In this particular case of EPDiff, denoted as EPDiff(Hl), one 
obtains the velocity u from the momentum m by an inversion 
of the elliptic Helmholtz operator (1 — a 2 A), with length scale 
a and Laplacian operator A. The EPDiff(Hl) equation with 
momentum definition (|5} therefore describes geodesic motion 
on the diffeomorphism group with respect to ||u||^-i, the 
norm of the fluid velocity. This velocity norm is recognized 
as being (twice) the "kinetic energy," when equation 
with momentum definition l|5} is interpreted as a model fluid 
equation, as in shallow water wave theory [27]-[29]. Although 
this operator does not appear in the class of operators used in 
CA (because it does not satisfy the p admissibility condition), 
this model exhibits a number of features which are highly 
relevant also in this case. In particular, singular momentum 
solutions emerge in the initial value problem for this model 
which behave as isolated waves, called solitons. In many ways, 
landmark matching in CA can be seen as generating a soliton 
dynamics between two sets of landmarks. 



B. Singular momentum solutions of EPDiff. 

In the 2D plane, EPDiff, Q, has weak singular momentum 
solutions that are expressed as [26], [28] 



N 

m(x, t) = ^ p a(t, s)5(x - Q (t, a)) ds . 

„_i J s 



(6) 



where s is a Lagrangian coordinate defined along a set of N 
curves in the plane moving with the flow by the equations 
x = Qa{t,s) and supported on the delta functions in the 
EPDiff solution (|6jl. Thus, the singular momentum solutions 
of EPDiff are vector valued curves supported on the delta 
functions in representing evolving "wavefronts" defined 
by the Lagrange-to-Euler map for their momentum. These 
solutions have the exact same form as the landmark solutions 
obtained in the previous section (with the straightforward 
extension of matching 2N curves instead of 2). 

Substituting the defining relation u = G * m into the singu- 
lar momentum solution (|6} yields the velocity representation 
for the wavefronts, as another superposition of integrals, 



N 

U(X, t)=J2 P a(t, S)G(X, Q a (f , s)) ds 
„_i J s 



(7) 



In the example of EPDiff(Hl), the Green's function G for the 
second order Helmholtz operator in (0 relates the velocity 
to the momentum. In this case, the velocity in the singular 
solution Q has a discontinuity in its first derivative (its slope) 
across each curve parameterized by s moving with the flow. 
Being discontinuities in the gradient of velocity that move 
along with the flow, these singular solutions for the velocity 
are classified as contact discontinuities in fluid mechanics. 
These contact discontinuities do not occur if the Green's kernel 
is sufficiently smooth (eg. a Gaussian kernel), as is typically 
used in CA. 

C. Lagrangian representation of the singular solutions of 
EPDiff. 

Substituting the singular momentum solution formula l|6} 
for s £ M 1 and its corresponding velocity Q into EPDiff 
(0, then integrating against a smooth test function implies the 
following Lagrangian wavefront equations 



d_ 

dt 



Q a (s,t) = fPb(s , ,t)G(Q a (s,t),Qb(s' 1 t))d s ' 

6=1 ^ 

r) N r 

^P Q (M) = E J *P.(M)P&(«',*) 



(8) 



d 



dQ a {s,t) 



G(Q a (s,t),Qb(s',t))ds' . 
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Thus, the momentum solution formula yields a closed set 
of integral partial differential equations given by (jSJl for the 
vector parameters Q a (s, t) and P (s, t) with a = 1,2 ... N. 
The dynamics (|8} for these parameters is canonically Hamil- 
tonian and geodesic in phase space. 

D. Relation between contact solutions of EPDiff(Hl) and 
solitons. 

As we have discussed, the weak solutions of EPDiff(Hl) 
represent the third of the three known types of fluid singular- 
ities: shocks, vortices and contacts. The key feature of these 
contacts is that they carry momentum; so the wavefront inter- 
actions they represent are collisions, in which momentum is 
exchanged. This is very reminiscent of the soliton paradigm in 
ID. And, indeed, in ID the singular solutions Q of EPDiff are 
true solitons that undergo elastic collisions and are solvable by 
the inverse scattering transform for an isospectral eigenvalue 
problem [9]. Besides describing wavefronts, this interaction of 
contacts applies in a variety of fluid situations ranging from 
solitons [9], to turbulence [12], [15]. The nonlocal elliptic 
solve u = G * m relating the momentum density m to the 
velocity u in the EPDiff(Hl) equation also appears in the 
theory of fully nonlinear shallow water waves [9], [21], [28], 
[29], [49]. 

The physical concept of momentum exchange is well un- 
derstood for nonlinear collisions of shallow water waves, 
especially in ID. Momentum exchange for EPDiff in ID is 
exhaustively studied in [16]. The corresponding momentum 
exchange processes (wavefront collisions) for EPDiff(Hl) in 
2D and 3D are studied in [28], [29]. These wavefront collisions 
show an interesting phenomenon. Namely, wavefront solutions 
of EPDiff(Hl) in 2D and 3D for which a faster wavefront 
obliquely overtakes a slower one result in the faster wavefront 
accelerating the slower one and reconnecting with it. Such 
wavefront reconnections are observed in nature. For example, 
observations from the Space Shuttle show internal wavefront 
reconnections occurring in the South China Sea [32]. This 
wave front reconnection phenomenon is illustrated in Figure 
[0 The key mathematical feature responsible for wave front 
reconnection is the nonlocal nonlinearity appearing in EPDiff. 

E. norm versus smooth kernels 

As noted before, the Green's kernels for the operators L 
used in CA are typically smoother than the inverse of the 
elliptic Helmholtz operator (1 — a 2 A) which corresponds to 
the i/^-model. A consequence of this is that the (variational) 
matching problems are always well-posed in CA, and their 
solutions are computationally feasible. 



To illustrate the differences introduced by using smooth 
Green's kernels, consider the case of a symmetric head-on 
collision of two particles (or landmarks) in ID. Under the H\ 
model, they will meet in finite time, then bounce back after 
exchanging their momenta. This is impossible with a smooth 
kernel, since the landmarks are carried by a diffeomorphic 
motion, and therefore cannot meet if they started from different 
positions. 

F. Metamorphoses 

Interestingly enough, the crossover behavior can be recov- 
ered by using metamorphoses, replacing A by A a in (|4}, 
because, for this model, particles are slightly disconnected 
from the diffeomorphic motion, and in this precise situation 
slightly ahead of it, allowing them to cross without creating 
a singularity. After crossing, the landmarks carry the dif- 
feomorphism the other way, letting the motion appear like 
a compression, then a bouncing back after the crossover. 
To illustrate this, a symmetric frontal shock between two 
"landmarks" has been simulated in the cases a ~ vs. a > 
the result is in figure [2] 

In 2D, the same behavior can be observed. With a smooth 
kernel, two expanding circles which have no intersection at 
time t = will have no intersection at all times. This is shown 
in figure [5] When the parameter a 2 is introduced, the circles 
intersect in finite time, the value of a 2 influencing their shapes 
before and after the collision (figure |4] and [5}. 

IV. Applications of EPDiff in CA. 

So far, CA has been related to the issue of comparing 
two geometric objects, and thus more concerned with the 
variational "boundary value" problem 0. However, as studied 
in [34], the initial value problem, associated to the integration 
of EPDiff, turns out to have very important consequences for 
applications. 

The main consequence of the Euler-Poincare analysis is 
that, when matching two geometric structures, the momentum 
at time t = contains all the required information for 
reconstructing the target from the template. This momen- 
tum therefore provides a template-centered coordinate system, 
which essentially encodes all possible deformations which can 
be applied to it. 

There is another important feature which has been observed 
in the landmark matching problem. Although the momentum 
is a priori of functional nature, as a result of the application 
of the operator L to the velocity field u t , in the landmark 
case, we found that it was characterized by a collection of K 
vectors in space, for a matching problem with K landmarks. 
Thus, the momentum has exactly the same dimension as 



Fig. 1 

A SINGLE COLLISION IS SHOWN TO PRODUCE RECONNECTION AS THE FASTER WAVE FRONT SEGMENT INITIALLY MOVING SOUTHWEST ALONG THE 
DIAGONAL EXPANDS, CURVES AND OBLIQUELY OVERTAKES THE SLOWER ONE, WHICH WAS INITIALLY MOVING RIGHTWARD (EAST). THIS 
RECONNECTION ILLUSTRATES ONE OF THE COLLISION RULES FOR SINGULAR SOLUTIONS OF THE TWO-DIMENSIONAL EPDlFF FLOW. SEE [29] FOR A 

COMPLETE TREATMENT. 



the matched structures, and there is no redundancy of the 
representation. This is generic in the sense that the final 
dimension of the momentum exactly adapts to the nature of the 
matching problem. For example, as stated in [34], [35], when 
matching two smooth scalar images, the momentum has to be 
normal to the level sets of the image, and therefore is uniquely 
characterized by its intensity. The momentum can therefore be 
seen as an algebraic image, which has the same dimension as 
the original image. This is because template matching brings 
an additional reduction to the original analysis of finding 
optimal paths between diffeomorphisms (which leads to the 
EPDiff equation). This reduction is because the solution must 
be modded out by the diffeomorphisms that leave the template 
invariant, which constrains the optimal solution. As a result, 
the momentum is a (locally) one-to-one representation of the 
targets, in this template-based coordinate system. It is therefore 
a non-redundant tool for representing deformations of the 
template. 

Besides being one-to-one, the other advantage of the mo- 
mentum representation is that it is linear in nature, being dual 
to the velocity vectors. Thus, linear combinations of either 



velocity fields, or momenta are meaningful mathematically 
and physically, provided they are applied to the same tem- 
plate. Thus the average of a collection of momenta, of their 
principal components, or time derivatives of momenta at a 
fixed template are all well-defined quantities. 

Finally, because they provide an efficient tool for repre- 
senting deformable data, momenta are perfectly suitable for 
modeling deformations. Any statistical model on momenta 
provides, after the integration of EPDiff, a statistical model 
on deformations, the advantage being that it is much easier 
to build, sample and estimate statistical models on a linear 
space. An illustration of this is provided in figure [6] in which 
the momentum was generated as positive uncorrected noise 
on the boundary of a disc, and EPDiff was integrated with 
this initial conditions. The Green's function we used for this 
is a Gaussian kernel, G(x,y) = exp(— \x — y\ 2 ). Figure 
shows the initial momentum of the geodesic path between 
two 3D sets of landmarks placed on two hippocampi. We also 
refer to the statistical experiments on PCA in momentum space 
presented in [58]. 



Fig. 2 

Deformation resulting from a head-on collision of two Gaussian landmarks. First row: y coordinates of the evolving 

LANDMARKS PLOTTED AGAINST TIME FOR cr = (LEFT) AND cr > (RIGHT): CONTACT REQUIRES AN INFINITE TIME WHEN cr = 0, AND A CROSSOVER 
IS OBSERVED IN THE SECOND CASE. IN THE SECOND ROW, THE GRID DEFORMATION IS PLOTTED FOR THE ASSOCIATED 2D DEFORMATION, IN THE 

EXACT MATCHING CASE (cr = 0); THE GRID GETS SQUEEZED WHILE THE DISTANCE BETWEEN THE LANDMARKS REDUCES. IN THIS CASE THE 
DEFORMATION IS EXACTLY CARRIED BY THE LANDMARKS. THE THIRD ROW IS WITH a > (METAMORPHOSIS): IN THIS CASE, THE LANDMARKS 
TRAVEL SLIGHTLY AHEAD OF THE DEFORMATION, AND CAN CROSS WITHOUT CREATING A SINGULARITY. BEFORE THE CROSSOVER, THE GRID 

CONTRACTS, THEN THE LANDMARKS DIVERGE AND THE GRID EXPANDS. 




Fig. 3 

Motion of two expanding circles through EPDiff with a smooth kernel. This corresponds to the case a = (exact matching). The 

CIRCLES BECOME INFINITELY CLOSE, WITHOUT CROSSING. 
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Fig. 4 

Motion of two expanding circles through metamorphoses with small u' 2 > 0. The circles intersect in final time, like the 

LANDMARKS IN THE HEAD-ON COLLISION OF FIGURe|5] 




Fig. 5 

Same as in figure|4] with a larger a 2 > 0. 




Fig. 6 

Three deformations of a disc by EPDiff under random initial conditions. These figures visualize the evolution of a disc under 
the evolution of epdiff, for initial uncorrelated noise momentum on its boundary. this shows how random deformations in 
momenta superpose linearly to produce a diffeomorphic change in shape. 



V. Conclusions 

We have identified momentum as a key concept in the 
representation of image data for CA, and discovered im- 
portant analogies with soliton dynamics. Future work will 
explore further applications of EPDiff. In particular, it will 
be interesting to see whether the exchange of momentum in 
the interactions of multiple outlines will become as useful a 



concept in computational image analysis as it is in soliton 
interaction dynamics. 
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Fig. 7 

Landmark matching for hippocampi. Two sets of landmarks have been manually defined on the hippocampi in the first row, the 
first one being considered as the template. the second row shows the initial momentum superimposed with the template (left) 
and alone (right). most of the momenta is concentrated at the head and tail of the hippocampus. data taken from the 

Biomedical Informatics Research Network IWWW . MisikN . MtlTI 
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